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THE INTERACTION AND PENETRATION OF GASEOUS JETS 
IN SUPERSONIC FLOW 

by 


R. C. Orth, J. A. Schetz, and F. S. Billig 

Johns Hopkins University 
Applied Physics Laboratory 
Silver Spring, Maryland 


SUMMARY 


Experiments and analysis of the penetration of a gaseous jet 
from a discrete orifice into a supersonic stream and the associated 
interaction with the external flow are discussed. Included in the 
interaction is the definition of the shape of the stand-off shock in 
the main flow. A governing parameter in this analysis is the degree 
of underexpansion of the jet which has been shown to be directly re- 
lated to the ratio of jet static pressure to the pressure downstream 
of a strong disturbance in the main stream. For underexpanded jets 
the location of the Mach disk which is paramount to the definition 
of the penetration has been explicitly given. For matched pressure 
injection or for the region of flow downstream of the Mach disk, a 
solid body drag model has been revised and extended to describe the 
penetration into a supersonic stream. 

The study encompassed the effects of free stream Mach number, 
hole size and shape and the molecular weight, Mach number, and pres- 
sure of the injectant. An exemplary case of the design of the 
injectors for a typical supersonic combustor is given. 
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NOMENCLATURE 


speed of sound 

jet cross-section area 

separation height 

drag coefficient, D/hds 

drag, or combustion chamber diameter 

initial jet diameter 

equivalent jet diameter at sonic conditions 
fuel/air equivalence ratio 
acceleration of gravity = 32.17 ft/sec a 
width of jet 
mass 

mass flow rate 
Mach number 

number of injection ports 
static pressure 

average pressure over jet cross section 
effective back pressure = f P' 

L a 

pu s /2 = dynamic pressure 
radius of jet 
radius of curvature 
equivalent body nose radius 
distance along jet axis; "s = s/d 
time 

temperature 

axial velocity in jet 
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velocity in undisturbed air stream 
initial jet velocity 

downstream distance from the injector port 
abscissa of center of Mach disk 
arbitrary downstream location, Eq.(19) 
normal distance from the plate surface 
ordinate of Mach disk 
ordinate of outer jet boundary 
total displacement of Mach disk, Eq.(13) 
angle of jet incidence and initial value 
limiting angle in Eq. (12) 

angle of jet incidence at edge of separation layer 
displaced height of shock inflection point 
ratio of specific heats 
boundary layer thickness ahead of jet 
shock stand-off distance 
inlet kinetic energy efficiency 
( 1 + / ( 1 + ) 
density 


Subscripts 

a = undisturbed mainstream conditions 

j = initial jet conditions 

js = initial conditions for supersonic injection 

t = total conditions 

oo = free stream conditions 
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Superscripts 

( )' = conditions following normal shock in main stream 

( )" = conditions corresponding to sonic point in jet 
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INTRODUCTION 


In recent years there has been considerable interest in the flow- 
field produced by the angular injection of fluid (liquid or gas) 'through 
discrete openings in a wall into an otherwise uniform supersonic stream. 

Past work in this area has been aimed primarily at aerodynamic control 
systems; 1-6 consequently, major emphasis has been placed on the pressure 
field produced on the surface near the injection station. The development 
of scramjet engines, which require fuel injection into a supersonic airflow, 
has underscored the importance of understanding the mechanism of jet pene- 
tration so that the resulting jet trajectory and fuel-air mixing can be 
predicted. 

The main contribution of prior work on penetration is the presen- 
tation of reliable experimental data and correlations of various injector 
and free-stream parameters. Most often used are the injectant-to-free- 
stream pressure and momentum ratios. 7 11 These correlations are based on 
data measured at stations more than ten injector diameters downstream from 
the injector (x/d ^ > 10). Thus, they provide estimates of the gross per- 
formance of the injectors, but the results are insufficient for a clear 
understanding of the jet penetration in the near downstream region 
(x/d < 10). 

The analytical work has generally been rather crude, involving 
one or more of the following untenable assumptions: 1) the injected flow 

turns sharply and then follows the surface; 2) the shape of the turning 
jet is semicylindrical; 3) the flow in the jet is homentropic; 4) mixing 
is instantaneous at the injection station and dissipates the transverse 
jet momentum; 5) small disturbance theory can predict the pressure field; 

6) the jet acts as a rigid vertical obstruction with some arbitrary height; 

7) the total injectant penetrant is the height of the first Mach disk 

in an underexpanded jet; and 8) the transverse momentum of the jet remains 
constant. 


This report describes the analytical and experimental work done 
at APL in the area of jet penetration over the past three years. Por- 
tions of the work have been previously described in abridged forms 
(Ref. 12-19). 


ANALYSIS 


JET TRAJECTORY ANALYSIS 

The interaction of a secondary jet and subsequent dispersion 
into the primary fluid is assumed to be a two-stage process. In the 
penetration stage, the jet retains its identity while being accelerated 
and turned in the flow direction of the primary fluid. The second stage 
is considered to be a coaxial turbulent mixing process and is not dis- 
cussed herein. We shall retain the assumption that has been used in 
subsonic analyses 8 that the jet acts as some "solid" body that emerges 
from the port and is bent downstream by drag and distorted in cross-section 
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by pressure differences on the front and back faces and viscous shear. 

Two general cases are considered, viz. , the "on design" case in which 
the pressures of the jet and the main stream are nearly matched, and 
the underexpanded jet. In the latter case, the turning and distortion 
of the jet is preceded by a rapid area expansion and normal shock. The 
analysis for this case is identical with the former except that the 
initial conditions (density, velocity, cross-sectional area, and angle) 
correspond to the conditions following the normal shock, and the height 
to the normal shock or Mach disk is superposed to the penetration depth. 
The horseshoe -shaped cross section of the turning jet which has been 
observed experimentally 7 ’ 8 is approximated as an ellipse, and the growth 
of the width is taken from subsonic measurements. 8 The drag coefficient 
on an element of the jet "body" is taken as that for an infinite cylinder 
at the local angle of incidence to the stream. 20 

Consider the forces acting on an infinitesimal length ds of the 
jet cross section, as shown in Fig. 1'. First, balancing forces along 
the jet axis, 


m(du/dt) = -mg sin CC + APA (1) 

where AP is the change in average static pressure acting over the cross 
section in going from s to s + ds. Then, balancing forces on an element 
of mass m in a direction perpendicular to the axis, 

D + mg cos a = mu 2 /R (2) 

where u = u(s) is the local jet velocity and R = R(s) is the radius of 
curvature of the jet centerline. Now, an increment of time dt is related 
to the length and velocity by dt = ds/u, which can be substituted into 
Eq. (1). Furthermore, ds = RcKX , and D = Cjjq a (sin 2 a) hds , where h is the 
local width of the jet measured in the plane perpendicular to the plane 
in the sketch in Fig. 1. Generally, the forces due to drag and inertia 
are large compared to gravity forces, and the latter can be neglected. 
With this assumption and the foregoing substitutions, Eqs. (1) and (2) 
become 

pu(du/ds) = -dP/ds (3) 

dCt/ds = -C hq sin 2 Ct/pu 2 A (4) 

L) cL 


Since puA is the mass flow in the jet, which remains constant at the 
initial value, Eq. (4) can be put in the form 


da 

ds 


C_hp u 2 sin 2 a 
D a a 

2 P j U j a(A j a/A)(P j /p) 


(5) 


When u is eliminated from Eq. (4) in this way, it is not necessary to 
consider Eq. (3) to obtain the solution if h(s) , A(s) , p (s) , and 
Cp(a) are known. The last item can be determined by assuming that 
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Fig. 1 SCHEMATIC REPRESENTATION OF JET CENTERLINE TRAJECTORY 



Fig. 2 VALUES OF THE INTEGRAL IN EQUATION (11) FOR VARIOUS FREE 
STREAM MACH NUMBERS 
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each element of the jet acts as an element of an infinite cylinder 
aligned at the local angle to the flow. Experimental data are avail- 
able for this situation, 20 and the following simplified curve fits are 
taken as an adequate representation of these data: 

7/2 

Ck = 1.2+ (M sin a) ' 0 < M sin a < 1 (6a) 

D v a ' — a- v ' 

C_. = 1.06 + 1.14 (M sin Ct) ^ M sin a > 1 (6b) 

D a a — 

From subsonic experiments , 8 

h = 2.25 d. + 0.22 s (7) 

J 

Reference 8 also suggests that the shape of the cross section be taken as 
an ellipse with a ratio of major to minor axes of 5 to 1. With this and 
Eq. (7), the A(s) can be determined from 

A(s) = TTh 2 / 20 (8) 


Furthermore , 

h(A/A. s ) = [2.25 + 0.22 (s/d . ) ] 3 /2. 5TTd^ 
Finally, combining Eqs. (9) and (5), 


da 


f 3 

ds 2.5 tt p.u. p. 

J J J 


(9) 


( 10 ) 


where s = s/d., and Cp(tt) is obtained from Eq. (6a) or (6b) for the appro- 
priate range ^of M a sin a. The only remaining unspecified quantity is 
p(s), which in this analysis is assumed to be equal to the initial value 
either a.t the jet orifice or just downstream of the jet shock for under- 
expanded injection (see next section). Using this assumption, Eq. (10) 
can be integrated to give 


0.88 (2. 5tt) f a 0 da' 

( p a u a a / p j u j a ) J a C D (a')sin y (a') 


(0. 22 s + 2.25) 4 - 2.25 4 


( 11 ) 


The integral in Eq. (11) is evaluated by choosing a lower limit for the 
angle of incidence of interest, say 10° . One then integrates from this 
lower limit to 

a, = sin 1 (12) 
la 


using Eq. (6a) and then from a to CLq using Eq. (6b). Typical results are 
given in Fig. 2 for various values of M„. For M a < 0.2, it is sufficient 
to take = 1.2. 
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Model of the Flowfield for Angular Gaseous Jet Injection 


In Fig. 3 are shown schematic representations of typical flow- 
fields that have been observed for angular injection of gaseous jets 
into a supersonic free stream. There are four separate cases, and 
the occurrence of each depends upon the size of boundary layer and 
separation zone ahead of the injector and whether the injection pres- 
sure is matched or underexpanded. In general, a shock in the main- 
stream is caused by the obstruction of the primary flow by the emerging 
jet. If the pressure rise on the surface caused by this shock would 
be greater than about three to one, the boundary' layer will separate, 31 
and an oblique shock off the separated region will intersect the stronger 
mainstream shock. When the separation zone is small (Fig. 3a and 3b), 
the jet is immediately exposed to the momentum of the freestream and 
is displaced. The amount of displacement is, of course, dependent on 
flow conditions (e.g., PjUj 2 /p a u a 2 and M a ) . In any event, the pressure 
field on the periphery of the emerging jet is extremely complicated. 

Since the relation of Pj to this pressure field controls the internal 
structure of the jet flow, as well as the size of the separation zone 
ahead of ; the jet zone, some estimate of the pressure field must be 
obtained. 


Even if the pressure field around the jet were known in detail, 
there is no simple method known for predicting its effect upon the 
internal flow in the jet. Some insight into this general flow situa- 
tion can be obtained by studying the much simpler case of a jet entering 
a quiescent environment. 10 s11 The governing parameter for the jet 
structure in that case is the ratio of Pj to the "back pressure," 
i.e. , P|j for the quiescent medium. With this ratio and y, the results 
of Refs. 10 and 11 can be used to determine yQ and Mj , respectively 
(Figs. 3b and 3d). To apply these results to the case of a jet ex- 
hausting into a moving medium, we introduce here the concept of an 
"effective back pressure" P^ . If the emerging jet is viewed as a 
roughly cylindrical obstacle, the static pressure distribution around 
the periphery would be expected to vary from a maximum at the forward 
stagnation point down to a low value in the "base" region behind the 
jet; Pb is expected to be somewhere between these extremes. As a 
corollary to this argument, there will be a value of Pj, such that 
the structure in the jet is a minimum. The situation is defined as a 
"matched" pressure case; cases with Pj > P^ are "underexpanded". 

Studies at the University of Maryland 21 ,22 indicated that the 
average pressure field around a cylindrical protuberance on a flat 
plate in a supersonic stream was approximately equal to two thirds 
of the pitot pressure in the stream §P ta '. 21,22 This criterion is 
used in the present work to specify matched-pressure injection as that 
in which the static pressure of the secondary jet (Pj) equals the 
effective back pressure P^ = I'Pta'* The studies also showed that 
for the case of underexpanded injection the distance to the first 
Mach disk was only weakly dependent on the height of the separation 
zone. 21 The importance of this fact is that the Mach disk location, and 
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(c) THIN BOUNDARY LAYER, UNDEREXPANDED (d) THICK BOUNDARY LAYER, UNDEREXPANDED 


Fig. 3 SECONDARY JET STRUCTURE 



Fig. 4 PENETRATION FOR TWO FREE STREAM CONDITIONS WITH THE 
SAME RATIO, Pj0 j 2 /p a u a 2 
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thus, the value of qj at that location, can be determined with the 
knowledge of Pj/Pj, • The trajectory calculation is started from 
that point for all cases of underexpanded injection. 

Since the pressure rise associated with matched-pressure in- 
jection is sufficient to separate the boundary layer in supersonic 
flow trajectory, calculations for this case will have to include an 
estimate of the drag in the separated zone and the separation height 
b. With these estimates and the assumption that b/d^ = ¥ sin [ (CX ^ + 0.2)12], 
Eq. (11) can be solved for O 2 , the angle of jet incidence at the edge 
of the separation layer. The trajectory in the freestream can then be 
calculated using the jet initial conditions and O 2 as the initial angle. 

Effects of Initial Conditions on Jet Trajectory 

We are now in a position to explore the effects of some simple 
variations of initial conditions on the fuel jet trajectory. First, 
consider holding pjuj 2 /p a u a 2 constant while varying M a . The results 
of such a comparison are shown in Fig. 4 for M a = 4.0 and M a « 1. 

The interesting, if not surprising, result is that slightly better 
penetration is achieved for M a = 4. This behavior can be traced to 
Fig. 2, which shows a lower value of the integral in Eq. (11) for the 
range 90 < a < 18°. This results in a steeper slope of the curve in 
this range. The cause of all this is that Cp is lower at high normal 
Mach numbers (M a since) than for transonic or subsonic values, as evi- 
denced by Eq, (6). Finally, we may observe the effect of Ctg on pene- 
tration; results for a typical case are given in Fig. 5. 

EXPERIMENTAL PROCEDURE AND TEST RESULTS 

The objective of the experimental program was to study the 
effects of certain parameters on the jet-freestream interaction, i.e. , 
the jet structure and penetration in the near downstream region. The 
variables examined were the injectant-to-freestream pressure and momen- 
tum ratios, injectant and freestream Mach numbers, and the shape of a 
single injection port. To simulate the design restrictions imposed on 
a fuel injector by a fixed ER requirement, injectant mass flow was held 
constant during the parametric variations. 

Parametric Variation of Injection Pressure and Mach Number 


Experiments were conducted on a flat plate mounted in line with 
a Mach 2.72 free air jet (Fig. 6). Cold Eb? and N 2 were injected normal 
to the airstream from an axisymmetric sonic or supersonic nozzle located 
1.75 in. from the leading edge of the plate. The supersonic injectors 
(Mj s = 1.31, 1.50 and 1.67) had 10-deg. half-angles. Hydrogen concen- 
tration profiles were determined by gas -sampling with a seven-point, 
movable sampling rake. The samples were analyzed on a Beckman GC-2A 
Gas Chromatograph, and hydrogen concentration was determined by the 
peak height method. 23 To verify the accuracy of the analysis the sum 
of the partial pressures of hydrogen and air was compared with the 
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PjUj 
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Fig. 5 PENETRATION FOR A TYPICAL CASE WITH INJECTION AT TWO 
ANGLES TO THE FREE STREAM 


Fig. 6 FLAT PLATE INJECTOR AND MACH 2.72 FREE JET NOZZLE 
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VERTICAL DISTANCE FROM 
NOZZLE EXIT PLANE TO MACH DISK (y 0 /dj*) 


Fig. 7 RATIO OF TOTAL INJECTANT PRESSURE TO FREE STREAM 
EFFECTIVE BACK PRESSURE VERSUS MACH DISK DISTANCE 


PROBE STATIONS 



DOWNSTREAM DISTANCE FROM INJECTOR CENTERLINE (x/dj ) 


Fig. 8 JET TRAJECTORY AND INJECTANT CONCENTRATION PROFILES 
FOR MATCHED-PRESSURE SONIC INJECTION 
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measured total pressure of the sample before analysis. The rms 
deviation of the summed partial pressures from the measured total 
pressure was 0.66% and the maximum deviation was 5.3%. This com- 
parison insured that there were no gross errors in the chromatograph 
analysis. Frequent calibrations using pure known gases were made to 
detect any contamination of the gas separation column and/or shifting 
of the sensitivity of the thermal conductivity cell which could cause 
compensating errors in the measured concentrations of the constitutent 
gases in an unknown mixture. On the basis of the foregoing proced- 
ures we believe that the actual value of hydrogen mole fraction of 
any one sample was within + 2% of the measured value. Pressure 
measurements were taken periodically in the air and injectant plenum 
chambers, at the exit wall of the free jet, and in the freestream by 
using the sampling probe as a pitot probe. 

Schlieren photographs were taken on ASA 3000 Polaroid film 
in a focal-plane camera. They were useful for analysis but, unfor- 
tunately, were not of good enough quality to present for publication. 
For matched-pressure sonic injection (Pj = P^) , there was no visual 

evidence of disturbance within the jet, which simply turned down- 
stream with an outer boundary very close to the limit of the M < 1 
flow region that is shown by the absence of shocks in front of the 
probes. Strong normal shocks (Mach disks, similar to Fig. 3) were 
present in both the underexpanded sonic jets and the supersonic jets. 

For underexpanded injection, the perpendicular distance from 
the flat plate to the Mach disk was scaled from the schlieren photo- 
graphs and plotted vs the ratio of jet total pressure to effective 
back pressure (P t j /P]-,) (Fig . 7). Crist et al. made a thorough inves- 
tigation of underexpanded jets in a quiescent medium . 10 They showed 
a relationship between the normalized distance to the first Mach disk 
(h/dj) and P t j/P a . (For a quiescent medium, p ta = p a or P oo.) Schetz, 
Hawkins, and Lehman 21 have extended this concept to injection into 
supersonic flow by plotting the vertical distance to the Mach disk, 
yg/dj , and the ratio yp/xg vs Pj /P a ' . Note that the correlating 
effective back pressure P^ used in Ref. 10 is P a ' rather than -fp^' 
and that the quantity gP a '/P ta ' varies from 1.18 to 1.33 for y = 1.4 
and 1.9 < M a < 4.5, which is the Mach number range of the data shown 
in Fig. 7. If the data were normalized on the P a ' basis they would 
be lowered with respect to the correlating line for the quiescent 
medium data, thus giving better correlation with some supersonic 
points but worse with others. Since there was no clear advantage in 
using P a ' as the correlating parameter, we chose to use §-P ta ' in our 
presentation to be consistent with the concept of the average pressure 
around a cylindrical obstruction. 

The structure of the emergent supersonic jet was similar to 
that of the underexpanded sonic jet. The supersonic injectors were 
designed to give the same static pressure at discharge as the sonic 
injector for a given mass flow; thus for all supersonic injectors 
Pj = P 5 . With coaxial supersonic streams the pressure balance does 
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occur when the static pressures are matched, and any discrepancies 
in flow direction from the present conical (rather than contoured) 
secondary jet are adjusted primarily through a series of weak oblique 
waves. However, with transverse injection the pressure balance occurs . 
when the static pressures behind strong disturbances in both streams 
are approximately matched. Whether the secondary jet is conical or con- 
toured is probably relatively insignificant. On this basis the super- 
sonic injectors tested were actually underex panded, and, therefore, 
the emergent jets expanded into the main flow until the Mach number 
was such as to produce the effective back pressure aft of the Mach 
disk. The values of yo/dj* for the supersonic cases were from 9 to 
30% higher than those for underexpanded sonic injection with an equal 
injectant stagnation pressure. There was also less downstream dis- 
placement of the Mach disk than with the equivalent underexpanded 
sonic jet; however, additional data, particularly at higher values 
of M j , is required to produce a correlation of xq, yo> and Ftj/^b f° r 
supersonic injection. 

The injectant concentration profiles (Fig. 8-10) were developed 
from a best fit to the experimental data of the types of concentration 
profiles reported by Zukoski and Spaid. 7 

Injector Shape Variations 


Two different experimental set-ups were used to study the effect 
of individual injector shapes on penetration. The first series of tests 
was conducted in a two-dimensional Mach 2.1 wind tunnel. Helium and 
air were injected transversely from a flat plate, using three different 
sonic injector configurations: a circular -ended, 1/4-in. by 1/16-in. 

slot with its long axis aligned with the main stream flow; the same 
slot with the long axis perpendicular to the main stream flow; and a 
0. 136-in. -diameter circular nozzle. The cross-sectional area of the 
circular nozzle was equal to that of the slot. Injectant-to-frees tream 
pressure ratio, injectant mass flow, and freestream Mach number were 
held constant, and penetration was determined from schlieren photo- 
graphs, using air as the injectant, and from gas samples using He as 
the injectant. 

The second series of tests was designed to study the jet 
trajectory in detail and to provide information on the shape of the 
jet cross-section in off-axis planes. The experiments were run in a 
Mach 2.7 free-jet with injection of hydrogen from the flat plate model 
shown in Fig. 11. The injection port diameter was increased from 0.14 
in. to 0.20 in. to aid in the visualization of the flowfield with 
schlieren. The elongated slot injector was geometrically similar to 
its previously tested counterpart. The important jet penetration 
parameters, p t j / and qj/q a , were closely matched in the two tests. 

At Mach 2.1 those ratios were 6.5 and 3.5, and at Mach 2.7 they were 
5.92 and 2.8. The sampling station was at x/d. = 9.2 in the former 
and at x/dj = 5 and 10 in the latter. Gas samples were taken at off- 
axis locations with a seven-point probe that could be translated in 
three directions. Improvements in gas chromatography techniques and 
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(a) 


Pj /P b = 1-48 qj /q, = 1.29 
P t j/P b = 2.80 dj =0.144 IN. 


PROBE STATIONS 



DOWNSTREAM DISTANCE FROM INJECTOR CENTERLINE (x/dj ) 
Pj/P 5 =2.73 dj /q a = 2.38 


P t j/P b = 5.18 dj = 0.112 IN. 
PROBE STATIONS 



0 1 2 3 4 5 6 7 8 9 10 11 

DOWNSTREAM DISTANCE FROM INJECTOR CENTERLINE (x/d,) 


Fig. 9 JET TRAJECTORY AND INJECTANT CONCENTRATION PROFILES 
FOR UNDER-EXPANDED SONIC INJECTION 
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DISTANCE FROM INJECTION PLANE (y/d *) DISTANCE FROM INJECTION PLANE (y/d | *) DISTANCE FROM INJECTION PLANE (y/'dp 



DOWNSTREAM DISTANCE FROM INJECTOR CENTERLINE (x/dp 



PROBE STATIONS 



DOWNSTREAM DISTANCE FROM INJEQTOR CENTERLINE (x/dp 


CTORY AND INJECT ANT CONCENTRATION 
FOR SUPERSONIC INJECTION 
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INJECTOR 
PROBE RAKE 


MACH 2.72 
NOZZLE 


Fig. 11 FLAT PLATE INJECTOR WITH THREE DIRECTIONAL MOVEABLE PROBE 
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equipment made it possible to attain a reproducibility of + 1% in 
hydrogen concentration by use of the integrated-signal analysis. 33 
This excellent repeatability is attributed partly to the precise 
matching of the flow in the two branches of the gas chromatograph, 
which was achieved by maintaining pressure regulation on double 
sonic orifices in each circuit. Another factor in the improvement 
was the use of electrical integration of the output from the thermal 
conductivity cell rather than ratios of peak heights. 

Figure 12 shows schlieren photographs of the penetration 
from the three injector configurations. The size of boundary layer 
separation increased as the disturbance in the main stream caused 
by the frontal area of the secondary jet increased. With the long 
axis of the rectangular nozzle aligned with the mainstream flow 
(Fig. 12a) , the ratio of the separation height to the diameter of 
the circular injector (b/dj) was 0.35. The values of b/dj for the 
circular nozzle (Fig. 12b) and for the rectangular nozzle (where 
dj is now the diameter of a circle of equivalent area perpendicular 
to the main stream. Fig. 12c) were 0.56 and 0.70, respectively. The 
jet structure was also altered by the shape of the injector, with the 
jet from the low-aspect -ratio injector (Fig. 12a) being the most de- 
formed. 


Figure 13 is a normalized scaled sketch of the fine jet 
structure deduced from each of the cases shown in Fig. 12. Although 
the shape of the injector strongly affects the shape of the fine struc- 
ture of the underexpanded secondary jet (Figs. 12 and 13), it does not 
appear that the use of low-aspect-ratio injectors will improve pene- 
tration. As shown by Fig. 13, the location of the normal shock in the 
jet (called the "Mach disk" for the circular case) and the penetration 
in the near -downstream region were practically the same for the three 
injector shapes tested. 

Figure 14 shows the measured hydrogen concentration contours 
at a downstream location of x/dj = 5. Except for minor variations 
in the contours there is very little difference in the jet penetration 
and lateral spreading in the near downstream region. 

The parametric variation of injectant pressure and Mach number 
established that the normal distance to the center of the Mach disk, 
or normal shock, in the secondary jet was a function of the jet-to-free- 
stream pressure ratio, which specified the degrees cf underexpansion 
of the secondary jet. The initial drag on the secondary jet, which 
will be determined by the injector shape, will affect the downstream 
displacement of the Mach disk and the ratio yg/xg. The experiments 
substantiate the argument that the injector shape will affect the shape 
of the secondary jet structure and the orientation of the Mach disk, 
but the normal location of the Mach disk, and thus the penetration in 
the near downstream region are not greatly influenced by the injector 
shape. Furthermore, since the loss of jet momentum PjUj S across the 
normal' shock is great, the orientation bf the Mach disk and the- direction 
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of the jet as it crosses the shock also have little effect on the addi- 
tional downstream penetration of the jet as a whole. 

COMPARISON OF THEORY WITH EXPERIMENTAL RESULTS 

To verify the trajectory theory, trajectories calculated by using 
Eq. (11) were compared with the experimentally determined concentration 
curves. For the case of matched pressure injection, qj was evaluated 
at the point of injection and Ctj = 90 . The integral was evaluated 
using Fig. 2 (it can also be done numerically) for 15° increments, and 
the trajectory was determined by evaluating the As" from increment to 
increment until the full range of 90° was covered. 

Because there was a relatively large separated boundary layer, 
the accuracy of the calculated trajectory rests heavily on the knowledge 
of the drag in that region. Although the details of the separated 
region were not precisely defined, it was assumed that the drag on the 
secondary jet can be approximated by the drag in the separated boundary 
layer of supersonic flow around a solid cylinder on a flat plate, and an 
estimate of Cp for the jet was calculated from the cylindrical surface 
pressure data presented in Ref. 22. The assumption that the drag on 
the fluid jet was equal to the drag on a solid cylinder resulted in a 
trajectory that was slightly below the experimentally determined maximum 
concentration locations (Fig. 8). However, in the absence of more defi- 
nitive data on the interaction of the jet and the separated flow, it is 
suitable for a conservative estimate of the jet trajectory. In the more 
practical cases of jet penetration, such as the design of a fuel injector, 
it may not be possible to measure the separation height. Therefore, a 
calculation that estimates the effect of the separated boundary layer 
on the trajectory of the jet but is based only on freestream conditions 
is desirable. To provide such a calculation, we compared C D for the jet 
as it penetrated a distance equal to b in a stream with no boundary 
layer, Cp = f (M a , a) as in the foregoing, and C D for the separated 
zone. Since the drag for the jet in the freestream is larger than the 
drag in the separated boundary layer, the penetrations in the two cases 
can be made equal if the freestream trajectory is vertically displaced 
such that the freestream drag times the vertical distance over which it 
acts is equal to the separation zone drag times the separation height. 

For the experimental data shown in Fig. 8, the offset distance (for 
which Cj} = 0) that was required to make the completely inviscid tra- 
jectory coincide with the trajectory shown was twice the compressible 
flow boundary -layer thickness as calculated by Pai 24 and Eckert and 
Drake. 25 Thus, in this procedure the trajectory for injection into 
a completely inviscid freestream is calculated, then it is displaced 
upward a vertical distance equal to twice the compressible flow boundary 
layer thickness. 

To calculate the trajectory of an underexpanded jet, the loca- 
tion of the Mach disk (xq , yQ) was determined from the experimental 
correlations. The normal displacement yg was taken directly from Fig. 7 
and the ratio of normal to downstream displacement yg/xg = ^ > was use d 
in accordance with the data from Refs. 13 and 21 shown in Fig. 15. The 
total displacement z , was calculated from the assumption 
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With this Mach number, q on the downstream side of the Mach disk was 
calculated from the normal shock equations, then the jet trajectory 
downstream of the Mach disk was calculated from Eq . (11). The proce - 
dure for calculating the tra jectories of the supersonic jets was the 
same, except that the actual measured values of xq and yQ were used, and 
the length of the supersonic portion of the nozzle was included in Yq. 

The trajectories for underexpanded jets calculated on this basis 
are characterized by a rapid penetration to a point where the Mach disk 
occurs and sharply reduces q, followed by a gradual increase in pene- 
tration that occurs downstream of the Mach disk is due to the low value 
of q/q a in this region. 


The location of the Mach disk and the resulting loss of energy 
is a strong function of P t j/P^. The method of calculating the trajectory 
for underexpanded jets and supersonic jets downstream from the Mach disk 
rests heavily on the previously described extension of the correlation 
of z vs Mj (Eq . 14) in the quiescent medium to the case having exter- 
nal flow. Instream pitot and static pressure measurements in the region 
of the Mack disk are needed to ascertain the validity of this method. 

Comparison of Some Penetration Correlations 


Two types of expressions for jet penetration were found, one 
based on the trajectory of maximum injectant concentration, and the other 
based on the trace of the outer boundary of the jet. Equation (11) and 
the following equations (15-17) are of the former type, whereas Eq . (18) 
is of the latter type. Figures 16 and 17 compare these correlations with 
data from two of our tests with underexpanded hydrogen injection into a 
Mach 2.72 airstream. Thus 


y/ d . 
J 

= [(q./q ) x/d ] 0,394 

J d J 

(Ref. 8) 

(15) 

y/d . 
J 

, , .0.434, ,, .0.333 

= (q./q ) (x/d ) 

J a J 

(Ref. 8) 

(16) 

y/ d . 
J 

, ,, r „ / » \ s ,0.533, ,-0.259 

= 1.64 [q .M /q M.)cos s a.J (x/d.) 

J a a j J J 

(Ref. 9) 

(17) 

y/dj 

- 1.68 [q .M /q M. ) cos^o; . J (x/d.) 

j a a J J J 

(Ref. 9) 

(18) 
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Equations (15) and (16), which were developed for a subsonic 
main flow, 8 agree reasonably well with the experimental values for 
x/ d j <8 but lead to unreasonably large penetrations farther downstream. 

The values of the exponent on (qj/qa) in these formulas are satisfactory 
for supersonic flow (shown by the correspondence of the trajectories at 
x/dj = 1.0), but the values of the exponent on (x/dj) are too high; they 
predict the slowly increasing penetration typical of subsonic flow rather 
than the rapid penetration, quick turning, and flat trajectory of the 
secondary jet in a supersonic freestream. 

Equation (17) also was base£ on data far downstream, (x/d* = 7.5 to 
72). The negative exponent on (x/dj) attempts to account for the unsymmetric 
mixing that occurs in this far downstream region and causes the maximum 
concentration to shift toward the flat plate. The correlation is unreal- 
istic for x/d^ < 6. 

Only Eq. (11) considers the degree of underexpansion as well as 
q i /q a . Of the four correlations it appears to describe best the location 
of the maximum concentration trajectory. Furthermore, because the struc- 
ture of the emergent jet is considered, it is the most realistic. 

The location of the jet boundary or maximum penetration is partly 
a matter of definition. Vranos and Nolan 0 defined the jet boundary as 
the y/dj value at which the mole fraction of the injectant is 0.005, and 
they used this definition in the development of Eq. (18). Another defi- 
nition of the boundary is the y/dj value at which the mole fraction has 
fallen to 10-20% of the peak value at the same x/dj in the plane of the 
jet centerline. The latter definition was more convenient for determining 
the jet boundary from the experimental concentration curves in the present 
work (Fig. 8-10). Compatibility between the two definitions is shown by 
the agreement between the jet boundary as calculated with Eq. (18) and 
the experimental values based on the latter criterion for H 2 injection 
into a Mach 2.7 2 air stream (Fig. 17). 

Since the foregoing comparisons have been concerned with the 
relative locations of the maximum concentration (y/dj) and the outer 
jet boundary ( y/dj), it is now important to examine the effect of 
q • / q a and the degree of underexpansion on the absolute penetration. 

Tne importance of evaluating injectors on this basis is that when the 
mass flow of the injector is specified, as by a required fuel equiv- 
alence ratio for example, the better injector will be the one that gives 
the greatest fuel penetration. As shown in Table 1 the effect of under- 
expansion, or increased injectant pressure, on the absolute penetration 
is small. When the injection pressure is increased by a factor of 2.73 
relative to the matched-pressure case, the increase in absolute pene- 
tration based on maximum concentration (y) is only 7 or 8%. This small 
increase in the absolute penetration of the maximum concentration is 
due to the unsymmetrical mixing of the injectant and the freestream 
which occurs closer to the injector for higher degrees of underexpansion. 

For fuel injectors this phenomenon may indeed be beneficial, because 
the fuel would be coming into more intimate contact with the air, 
resulting in more rapid ignition. Since this is due to co-current 
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Table 1 


Effect of Injection Pressure on Absolute 
Penetration Measured at x = 1.0 In. 


V p b 

d , in. 

x/d . 
1 

y, in. 

y, in. 

0.75 

0.204 

4.90 

0.286 

0.530 

1.48 

0.144 

6.94 

0.288 

0.518 

2.73 

0.112 

8.92 

0.308 

0.588 
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mixing, a quantitative evaluation of this effect is beyond the scope of 
this work. Comparison of the outer jet boundaries (y ) in Table 1 
shows a 10% increase in penetration for Pj/Pb = 2.73. For either y 
or y , an increase in Pj/P^ from 0.75 to 1.48 had no appreciable effect. 
The general conclusion from Table 1 is that placing a large demand on 
a fuel pressurization or pumping system to increase fuel pressure for 
sonic injection will have only a small payoff. 

SUPERSONIC INJECTION 

One clear result of the analytical and experimental studies 
presented here and also from previously reported experimental studies 
is the obvious importance of qj/q a on penetration. Although increasing 
P t j increases this ratio, the effects are not large for P t . below the 
pressure match and are quite small for underexpanded injection. Accepting 
the result that Pj should be the matched pressure so that a relatively 
shock-free jet flow occurs, then it is instructive to explore another 
possibility for increasing q^/q a and hence penetration, viz., varying 
Mj . Since q = pyif /2, increasing Mj has the desired effect of increasing 
qj for the same static pressure level. 

Equation (11) can be presented in the same form as Eqs. (15) -(17) 
as shown below: 

y/dj = [(qj /q a ) (x/dj)] 0,4 (19) 

Consider the penetration of a given amount of fluid per unit time by nor- 
mal, sonic injection into a supersonic airstream. One can write 
q j = Pj u* s /2 = n>Uj/2Ai. Take as the standard for comparison a case 
with Pj = Pb- As opposed to this, take a supersonic case, denoted by 
subscript s, with an injection pressure Pj = P^. Now using Eq. (19), 
y s /y, at a given distance x downstream of the injection station, can be 
determined as 


(y/d.) (q, /q ) (x/d ) °' 4 

(y/df) = (q'j/q ) C?/d*) ( 19a : 

J J d J 

With x = x and q = q , this becomes 
s a as 

* 0.6 0.4 

y s /y = (d js /d j> (20) 
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We can now assess the magnitude of any advantage to be 
however, we must determine d.(M.) for a constant mass 
P . . Using ■*' ^ 
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flow and 


p = p* (1 + M s /5) 5/2 
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( 21 ) 


it can be shown that 
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A /A* = g°- V 1 
s 

where % = 1.2/(1 + 0.2M 3 ). Taking sonic injection as the standard case 

(d. /d’;) 0 ' 6 = §°* 15 M-°- 3 (23) 

js j' 


giving finally 


y s /y 


p 0. 15, ,0. 5 


(24) 


Figure 17 shows that experimental values of y s /y, the ratio of the peak 
concentration from supersonic injection to that for matched-pressure sonic 
injection, are in good agreement with the theoretical effect predicted 
for the particular case of matching of mass flow and static pressure at 
the injector exit. For Mj s = 1.67, the increase in penetration is 25% 

(ys/y = 1-25). Since an increase Hj S for a matched-pressure jet must be 
paid for by increasing P t j and since, as earlier stated, a more appropriate 
matching parameter of transverse supersonic streams is the static pressure 
behind strong disturbances in both streams, a comparison of the penetration 
of supersonic and sonic jets with equal P t j and mj was made (Fig. 18) . With 
these restrictions on the theoretical relationship, theoretical penetration 
gain peaks at Mj = (5)2. Unfortunately, the data scatter for this plot is 
large, and the variations in measured penetration ratios at different down- 
stream locations show no discernible pattern, so that further experimentation 
is needed to verify the theory. However, Figs. 17 and 18 do show a gain in 
penetration with supersonic injection. 


INTERACTION SHOCK SHAPE FOR TRANSVERSE INJECTION 


The effects of injector shape on jet penetration were further ex- 
plored by studying the interaction shock shapes produced by jets from diff- 
erently shaped orifices. Previous attempts by others (e.g.. Refs. 26 
through 28) to predict the shape of the interaction shock accompanying 
transverse injection across a uniform supersonic stream have been based 
on the concept of some "equivalent solid body," using results for shock 
shapes past blunt bodies. References 26 through 28 used the well-known 
blast-wave theory. The blast-wave theory is simple and convenient to 
apply, but suffers from the disadvantages that it cannot predict the shock 
stand-off distance and does not, in general, provide very good estimates 
of shock shapes past solid bodies. Recently, very simple expressions for 
shock shape and stand-off distance for blunt bodies have been shown to 
provide excellent predictions of both quantities (Ref. 29). 

The concept of an "equivalent solid body" is used in conjunction 
with simple expressions for the shock shape around blunt solid bodies 
to predict shock shapes produced by transverse injection into a super- 
sonic stream. To apply the shock shape expression (Ref. 29) it is 
necessary to determine some equivalent for the solid body nose radius 
in the transverse jet case. The one length scale that is directly known 
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is the radius of the injection port. The injection port radius by itself 
is not, however, a suitable choice for the "body" nose radius as it does 
not measure the size of the obstruction presented to the main flow. One 
can simply employ the equivalent nose radius given by the blast wave 
theory, i.e. , 




m . 

J 

Pa u a 


1. 

2 


(25) 


with the shock shape formula of Ref. 29. An alternative approach can be 
developed based on the work of Refs. (12), (21) and (30). Here the height 
of the Mach disk, h, can be chosen as a measure of the obstruction and 
used as the "body" radius. From the results of the cited references, a 
simple formula for this height may be written for cases with y^ = 1.4: 


h s R 


P . 
J 


-b 2 - P b 


(26) 


where K = 1.0 for Mj = 1.0. Both possible choices (Eq. (25) and Eq. (26)) 
will be tested by comparison with experiment. 


With these results, the formula of Ref. 29 applied to the present 
problem can be written; 
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where A/R is the nondimens ional stand-off distance given as 
B 

A/R = 0.143 exp [3.24/M 2 ], 

B a 


(28) 


and R/Rg is the nondimens ional radius of curvature of the shock at the nose 


R/R c = 1.143 exp [0.54/(M -l) 1 ' 2 ] 
■d a 


(29) 


Observations of schlieren photographs (Ref. 6 and Fig. 12) of trans- 
verse jet injection show that the interaction shock has an inflection point 
that is displaced above the solid surface. This can be interpreted to mean 
that the nose of the "equivalent solid body" is located above the surface. 

A good prediction of this displacement is necessary to provide an accurate 
picture of the shock shape in the near field. The displaced height of the 
inflection point is denoted as 8 , and it has been possible to correlate this 
height, when nondimensionalized, with the non-dimensional equivalent body 
radius, i.e. , 8 /r^ vs. Rg ^/r j , as can be seen in Fig. 19. 
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The predictions of this analysis have been compared to experimental 
results for the following three cases shown in Figs. 20 through 22, respec- 
tively: (a) sonic air injection through a round hole into a Mach 2.1 air 

stream; (b) sonic H 2 injection through a round hole into a Mach 2.7 air 
stream; and (c) sonic air injection with a 4 by 1 slot aligned with and 
across a Mach 2.1 air stream. Agreement is reasonably good in all cases, 
with the "body" radius of Eq. (26) providing a better prediction, especially 
for the hydrogen injection case. A result of significant practical importance 
emerges from this study. The analysis predicts that the shock shape is not 
a function of injector port shape; this is confirmed by experiment as seen 
in Fig. 22. It is important to note, however, that the experiments permitted 
only the planform of the interaction shock to be determined. Some difference 
in shock shape in planes other than that including the jet trajectory might 
exist, but would go undetected here. Nevertheless, the fact that the plan- 
form of the shock shape is insensitive to injector port shape is significant 
and adds further evidence to the conclusion that jet penetration is only 
weakly dependent on injector shape. 

APPLICATION TO SCRAMJET DESIGN 

At this point, it is informative to study the design problem of 
fuel injection in the combustion chamber of a hypersonic ramjet. Take a 
typical case as a q^ = 5000 psi trajectory with a velocity decrement from 
freestream to combustor inlet of 1500 fps and an inlet kinetic energy effi- 
ciency of 97%. Consider gaseous hydrogen with a total temperature of 294°K 
as the fuel, and take stoichiometric proportions of fuel and air. Let us 
investigate a flight Mach number of 6.0. For this case, an rigg = 0.97 trans- 
lates to a total pressure recovery of 51%. The 1500-fps velocity decrement 
gives a combustor Mach number of 2.2, and this results in a combustor stag- 
nation pressure of 1107 psia, a static pressure of 104 psia, q a = 351 psi, 
and (P ta ' ) = 695 psia. Estimating that "matched" pressure injection will 
occur, for sonic injection, at Pj = §(Pta') gives Pi = 464 psia, and 
q i = (y/2) (1) (464) = 325 psi. Thus, qj / q a = 325/351 = 0.925. 

For a circular combustion chamber of diameter D (ft), the air mass 
flow is 595 lb/sec. For stoichiometric hydrogen in air (0.0292 lb ILj/lb 
air) and n round injection holes of diameter dj , there results 

i 

(0. 0292) (595 D 2 ) =n(0.14)P t (TTd 2 /4) (144)T t 2 (30) 

j j 

or dj/D = 0.062 for n = 10. If we now select some axial distance downstream 
from the injection station, the fuel penetration can be displayed on a cross- 
section of the combustion chamber. Noting that for qj/<la ^ as i- n this 
case, the total fuel penetration is achieved in a short distance downstream 
(see Fig. 8), we select x = 4 dj as the axial station for display (open 
circles in Fig. 23a). (No attempt to show the distortion of the jet cross 
section has been made.) The striking result is that rather poor coverage 
of the air flow in the combustor has been achieved, and a low combustion 
efficiency would surely result. Therefore, one looks to see how this 
situation can be improved. Since the penetration is a function of the 
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Fig. 20 SHOCK SHAPE FOR SONIC AIR INJECTION 
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initial hole size and the hole size varies with the number of injectors, 
it would seem profitable to reduce the number of injectors. The results 
achieved by reducing the number of injectors from 10 to 4 are shown as 
open circles in Fig. 23b. Indeed, the fuel does penetrate somewhat farther 
into the air stream; however, larger air gaps are left around the periphery 
and one could say that the coverage of the air stream is, if anything, 
poorer in this case than with the 10-hole injector. 

It is instructive to consider briefly the effect of the choice of 
fuel on this question. Consider propane (CoHg) , which has a molecular 
weight of 44, a ratio of specific heats of 1.16, and stoichiometric fuel- 
air ratio of 0.064. The calculated dj/D is 0.044 for 10 holes and stoichio 
metric injection. Furthermore, the momentum ratio will be less than for 
hydrogen injection. Thus, since the initial diameter and the momentum 
ratio are both reduced, the penetration pattern for propane would be 
distinctly worse than that for hydrogen. The low density of hydrogen is, 
therefore, a real advantage with regard to penetration, since it keeps 
the initial jet diameter large. 

Returning now to hydrogen and noting the results of Table 1, it 
can be seen that some advantage can be gained by over-pressuring the jet. 

It does not, however, appear profitable to consider a pressure ratio 
greater than 2.0, since the gain in penetration diminishes rapidly, and 
P t j at an overpressure of 2.0 is already 1758 psia. The penetration 
pattern obtained in this way is shown as solid circles in Figs. 23a and 
23b. Although this is certainly an improvement, the 10-hole injector is 
still superior, and the pattern still leaves much to be desired. We saw 
in Figs. 18a and 18b that some gain in penetration could be obtained by 
supersonic injection but even these gains are not too significant. For 
example, if Pj s = Pj = 464 psia and mj s = mj but Mj is increased to 2, 

the penetration would increase by about 30% (Fig. 18a), but the required 
P t j s would increase from 878 psia to 3630 psia. On the other hand if we 
put the more reasonable constraint of P(-j s = P t j = 878 psia with m j s = mj 

and let Mj = 2 then the increase is but 15% (Fig. 18b) . Another factor 
that could have some effect on the fuel-air contact is the ratio of dj g /dj. 
For Pj s = Pj this ratio decreases with increasing Mj and is equal to 
0.639 at Mj =2 which probably is detrimental to mixing. On the other hand 
for Pm s = P t ^ the ratio d^ s /dj increases with increasing M^ and is equal 
to 1.30 at Mj = 2. 

One means of improving the over-all penetration pattern is to alter 
the shape of the combustion chamber itself. Figure 23c shows an elliptical 
(3:1) cross section of the same area as the circular cross section. Only 
the points corresponding to an overpressured jet are shown. Clearly, a 
much better fuel distribution is obtained in this way. (One disadvantage 
of this change is a 23% increase in the surface area of the combustor and 
thus an increase in heat transfer and skin-friction drag.) 

Now let us look at flight Mach number of 10 on the trajectory 
described previously. Here T^e = 97% gives a pressure recovery of 19%. 

We further get M a = 3.9 and qj/q a = 0.67. Matched pressure injection 
gives P t j = 750 psia. For stoichiometric proportions, 10 holes require 
dj/D = 0.042. If we envision the same combustor operating at both the 
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Fig. 23 COMBUSTION CHAMBERS FOR MACH 6 FLIGHT OF A 

SCRAMJET WITH STOICHIOMETRIC HYDROGEN INJECTED 



Mach 6 and Mach 10 flight conditions, then one could use only five of the 
holes at Mach 10. Since the diameter is then the same and qj/q a is very 
nearly the same for the Mach 10 case, the penetration pattern would be as 
shown in Fig. 23a or 23c using half of the holes. 

Cone lusions 


From the experimental studies it can be concluded that, as assumed 
in the trajectory analysis, the downstream region where the injectant 
behaves as a discrete jet with negligible mixing is quite short (x/ dj < 10). 

To describe accurately the jet behavior of the injectant farther downstream, 
coaxial mixing must be taken into account. 

The use of a "solid-body” drag model to describe the penetration 
of a fluid jet has been re-examined and extended to include the case of a 
supersonic external stream. A simple analysis was developed which relaxed 
most of the assumptions common to previous work. The most important step 
was the removal of the restrictive assumption that the normal component of 
the jet momentum is conserved. With this restriction removed, good agree- 
ment between the theory and experimental data was obtained using local 
values of the drag coefficient for an inclined cylinder. For underexpanded 
jets exhausting into the mainstream, the shock structure caused by under- 
expansion in the jet itself is the predominant factor governing penetra- 
tion. Introduction of an "effective back pressure" enabled the develop- 
ment of a Mach disk correlation analogous to the well-documented but much 
simpler case of underexpanded jets exhausting into a quiescent atmosphere. 

With this correlation the conditions of the jet just downstream of the Mach 
disk can be approximated and the resulting trajectories, which agreed closely 
with experimental results, are calculated. 

Several semiempirical penetration and trajectory correlations were 
compared with the experimental results obtained herein. They were based on 
other measurements taken far downstream (10 < x/d'j < 70). The trajectory 
calculation developed herein best described the location of the jet trajec- 
tory. With respect to the outer boundary of the fuel jet, the correlation 
by Vranos and Nolan agrees well with the near downstream data from the 
present studies. 

While the nondimensionalized penetration (y/dj) for sonic injection 
is governed by the degree of underexpansion and the initial ratio of jet-to- 
freestream momentum, variation in the absolute penetration due to unde rex pans ion 
is small if jets with equal mass flow are compared. The analytical analysis 
indicated that result and the experimentally measured concentration centerline 
profiles and crossectional contours verified the conclusion. In fact, due to 
the losses involved as the underexpanded jet crossed the normal shock (Mach 
disk), the absolute penetration for a jet with the underexpansion ratio of 
1.48 was not appreciably different (possibly even slightly less) than that 
of a jet whose injection pressure matched the effective back pressure of 
the supersonic freestream. A further increase to a 2.73 factor on under- 
expansion caused only an 8-10% gain in penetration. Thus, little is to be 
gained by increasing the injection pressure above the level of the effective 
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back pressure. Penetration can be increased by going to supersonic injec- 
tion. A theoretical comparison of penetration from a supersonic jet and 
a sonic jet of the same total iniectant pressure and mass flow showed a 
maximum penetration at Mj = (5)1/2. The experimentally measured penetra- 
tion at Mj = 1.67 was greater than that for Mj = 1.50 and in both cases it 
was 10 to 15% better than that for sonic injection. However, the scatter 
of the experimental points on this plot was large, and additional experi- 
mentation at higher values of Mj must be performed to verify that increas- 
ing Mj beyond (5)1/2 does not improve penetration. 

Absolute penetration is at best only weakly effected by injector 
hole shape. The momentum loss incurred by an underexpanded injectant as 
it crosses the Mach disk practically cancels any penetration differences 
caused by injector shape. Local concentration differences occur but these 
have little effect on the additional downstream penetration of the jet as 
a whole. 


The secondary effects of injector shape on penetration were further 
emphasized by the development of an expression for the interaction shock 
shape. The concept of an "equivalent solid body" with radius equal to 
the Mach disk height gave good results when used in conjunction with simple 
expressions for the shock shapes around blunt bodies. The agreement between 
the calculated shock shapes and those from schlieren photographs of the 
interaction shock caused by transverse injection was excellent. Finally, 
the analysis has been applied to the question of the design of a fuel injec- 
tion system for a supersonic combustion ramjet. Low qj/q a values are gen- 
erally to be expected for real flight cases. Under these circumstances, 
some underexpansion of the fuel jets is desirable, but it was shown that 
the attainment of a roughly uniform fuel air ratio throughout the cross 
section of a circular duct by fuel jet penetration is unlikely. 
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